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Microcanonical Thermodynamics 1] allows the application of Statistical Mechanics to finite and 
even small systems. As surface effects cannot be scaled away, one has to be careful with the standard 
arguments of splitting a system into two or bringing two systems into thermal contact with energy 
or particle exchange: Not only the volume part of the entropy must be considered. The addition 
of any other macroscopic constraint like a dividing surface, or the enforcement of gradients of the 
energy /particle reduce the entropy. As will be shown here, when removing such constraint in regions 
of a negative heat capacity, the system may even relax under a flow of heat against the temperature 
slope. Thus Clausius formulation of the Second Law: "Heat always flows from hot to cold" can be 
violated. However, the Second Law is still satisfied and the total Boltzmann-entropy is rising. In 
the final chapter the general microscopic mechanism leading to the convexity of the microcanonical 
entropy at phase separation is discussed. This is explained for the liquid-gas and the solid-liquid 
transition. 



I. INTRODUCTION 

In conventional (extensive) thermodynamics the ther- 
mal equilibrium of two systems is established by bring- 
ing them into thermal contact which allows free energy 
exchange 1 12|. Equilibrium is established when the total 
entropy 

Stotal = S\+ S2 (1) 

is maximal: 

dS to tai = dSi + dS 2 = 0. (2) 

Under an energy flux AE^i from 2 — > 1 the total en- 
tropy changes by 

AS total = [fix - /3 2 )A-E 2 ->i (3) 
(3 = dS/dE = I (4) 

Consequently, a maximum of Stotai{E = Ex +E2) will be 
approached when 

sign(AS tota/ ) = sign(T 2 - T 1 )sign(A£; 2 ^ 1 ) > (5) 

From here Clausius' first formulation of the Second Law 
follows: "Heat always flows from hot to cold". Essential 
for this conclusion is the additivity of S under the split 
(eqlQ . Temperature is an appropriate control parameter 
for extensive systems. 
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II. STABILITY AGAINST SPONTANEOUS 
ENERGY-GRADIENTS 

For a small or a very large self-gravitating system, c.f. 
PJ, additivity and extensivity of S, and also of E, is not 
given. In fact this is the main reason to develop this new 
and extended version of thermodynamics. At phase sepa- 
ration the microcanonical caloric curve T(E) is backbend- 

ing. Here the heat capacity C — — (ffj) 2 /fiJ7 becomes 
negative, the curvature §^ is positive, and S(E) is con- 
vex. This is the generic signal of a phase transition of 
first order [lj. The Potts-model illuminates in a most 
simple example the occurrence of a backbending caloric 
curve 0E3. A typical plot of s(e) = S(E = Ne)/N in 
the region of phase separation is shown in figQ. Sec- 
tion [ffl] discusses the general microscopic reasons for the 
convexity. 

Let us split the system into two pieces a + b by a divid- 
ing surface, with half the number of particles each. The 
dividing surface is purely geometrical. It exists only as 
long as the two pieces can be distinguished by their dif- 
ferent energy/particle e a and e&. If the constraint on the 
difference et — e a is fully relaxed and — e a can fluctuate 
freely at fixed e 2 = (e a + ef,)/2, the dividing surface is 
assumed to have no further physical effect on the system. 
Constraining the energy-difference e& — e a = Ae between 
the two, reduces the number of free, unconstrained de- 
grees of freedom and reduces the entropy by — AE. 

In the convex region (upwards concave like y — x 2 ) of 
s(e) this becomes especially evident: at e = e 2 ~ e m im 
e a = e 2 — Ae/2 and = e 2 + Ae/2, the simple algebraic 
sum of the individual entropies s sum = (s a + s&)/2 leads 
to s sum > s 2 . However, e S2 is already the sum over 
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FIG. 1: Potts model (q = 10) on a 50*50-lattice with periodic 
boundary conditions in the region of phase separation. At the 
energy e\ per lattice point the system is in the pure ordered 
phase, at eg in the pure disordered phase. At e a little above 
d the temperature T a = 1//3 is higher than T2 and even more 
than Tb at eb a little below ez . At e a the system separates into 
bubbles of disordered phase embedded in the ordered phase 
or at e& into droplets of ordered phase within the disordered 
one. If we combine two equal systems: one with the energy 
per lattice site e a — ei + Ae and at the temperature T a , the 
other with the energy e& = eg — Ae and at the temperature 
Tb < T a , and allowing for free energy exchange, then the sys- 
tems will equilibrize at energy ei with a rise of its entropy. 
The temperature T a drops(cooling) and energy (heat) flows 
(on average) from b —> a. I.e.: Heat flows from cold to hot! 
Clausius formulation of the Second Law is violated. This is 
well known for self-gravitating systems. However, this is not 
a peculiarity of only gravitating systems! It is the generic sit- 
uation at phase separations within classical thermodynamics 
even for systems with short-range coupling and has nothing 
to do with long range interactions. 



all possible configurations at e 2 including the split ones 
with e& — e a = Ae. Evidently, the enforced energetic 
split of the system into two pieces with prescribed specific 
energies eb > e a destroys configurations with different or 
vanishing — e a . This reduces the partition sum. It 
costs some additional "surface" or constraining entropy 
AE > to maintain the imbalance of the energy, such 



From this little exercise we learn that the microcanon- 
ical entropy per particle s(e,N,a) of a finite system de- 
pends not only on the energy per particle e = E/N 
and possibly still on the particle number N (for a non- 
extensive system) but also on the boundary (the con- 
tainer), or any other non-extensive external constraint, 
via £ = Na. If one splits a given configuration like in 
eq.Q one must not forget the change — AS in the surface 
entropy. 

The correct entropy balance, before and after estab- 
lishing the energetic split eb > e a of the system, is 



Safter ^befo 



Sb 



S2 — Act < 



(6) 



even though the sum of the first two terms is positive. 

In the inverse direction: Only by relaxing the con- 
straint and by allowing, on average, for an energy-flux 
(Aeft-nj > 0) opposite to T a — Tb > 0, against the 
temperature- gradient (slope), but in the direction of the 
energy-slope, an increase of stotai — * S2 is obtained. 
Then, on average both energies e a and eb become equal 
= ei and Act = 0. Or in other words, starting with an 
energy difference Ae — eb — e a > the allowance of un- 
constrained energy fluctuations around Ae = between 
the two subsystems implies the disappearance of Act and 
in general leads to an increase of the total entropy. This 
is consistent with the naive picture of an energy equilibra- 
tion. Thus Clausius' "Energy flows always from hot to 
cold" is violated. Of course this shows again that unlike 
to extensive thermodynamics the temperature is not an 
appropriate control parameter in non- extensive systems. 

In the thermodynamic limit N — > 00 of a system with 
short-range coupling AS ~ N 2/3 , AS/TV = Act oc N~ 1/3 
must go to due to van Hove's theorem. 



III. THE ORIGIN OF THE CONVEXITY OF 
S(E) AT PHASE SEPARATION 

Here I discuss the general microscopic mechanism lead- 
ing to the appearance of a convex intruder in S(E, V, ■ ■ ■ ) . 
This is the signal of phase transitions of first order and 
of phase-separation within the microcanonical ensemble. 

I assume the system is classical and obeys the Hamil- 
tonian: 



H = 



E 

i<j=l,JV 



Pi r 

2m 



+ V lnt (r, 



(7) 



that S 2 > S s 



AS. 



In this case the system is controlled by energy and vol- 
ume. 



A. Liquid— gas transition 

The microcanonical sum of states or partition sum is: 
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W(E, N, V) = *» u J *» ,0 - £ *L _ V**i? i} ] 

N\T(3N/2)(2irh 2 ) 3N / 2 J v n V n \ E - E J U 
= W ld ~ 9as (E - E ,N,V) x W mt (E - E Q ,N,V) 

= e [s"-'»+s««] (9) 

yJV E (3N-2)/2 yiN 3/2 

1 ' ' J N!T(3N/2) {2itH 2 ) 3N / 2 1 j 



V is the spatial volume. E'o = min y mt {r } is the energy 
of the ground-state of the system. 

The separation of W(E, N, V) into W ld - gas and W mt 
is the microcanonical analogue of the split of the canoni- 
cal partition sum into a kinetic part and a configuration 
part: 



JV 



Z ( T ) = T7T 



V 



M \2irh 2 



mT\ 3N/2 fd™ 



V N 



n v""{ 

e t 



(12) 



In the thermodynamic limit the order parameter of the 
(homogeneous) liquid-gas transition is the density. The 
transition is linked to a condensation of the system to- 
wards a larger density controlled by pressure. For a finite 
system we expect the analogue. However, here controlled 
by the constant available system volume V. If it is larger 
than the eigen-volume Vo of the condensate, the system 
does not fill the volume V at low E. N — 1 internal co- 
ordinates are limited to Vo. Only the center of mass of 
the droplet can move freely in V. The system does not 
fill the 3N-configuration space V N . Only a stripe of the 

1 /3 

width V in 3(N — 1) dimensions of the total 3-N dim 
space and with the length V 1 / 3 in the remaining 3 dimen- 
sions of center of mass motion is populated. The system 
is non-homogeneous even though it is equilibrized and, at 
low energies, internally in the single liquid phase. It is 
not characterized by an intensive homogeneous density. 
In fact, W int (E - E , N, V) can be written as: 



W int (E-E ,N,V) 



V(E) 



V 



jY 



N def 



d 3N r 



int 



[V{E)] 
[ E-E 



< 1 



(3JV-2)/2 



(13) 



S mt (E-E ,N,V) = N\n 



V 



(14) 
< (15) 



In eas. l|llll4|) the integration is only over regions in 3N- 
dim coordinate space where V mt {ri} is smaller than 
E, where the big brackets are positive. The volume 
V N (E) < V N is the accessible part of the spatial volume 
outside of the forbidden 3N-dim regions with V mt {ri} > 
E. I.e. the eigen-volume of the condensate (droplets) at 
the given energy. V N (E) has the limiting values: 



[V(E)] 



N 



V 



N 



Vo W_1 V for E 



for E in the gas phase 
E 



So that W mt (E ~E ,N, V) and S int (E -E ,N, V) have 
the property: 

W int {E-E ) < 1, S mt {E-E Q ,N) < 

{ 1 E > V mt 



E — > Eq 



S in \E - E Q ) 



l n {[to.]N-i} <0 E^E Q 



(16) 



(17) 



All physical details are encrypted in W mt (E — 
E Q ,N,V) or S int (E - Eq,N,V) alias iVln[V(^)], c.f. 
eas. l(T3HT7|) : If the energy is high the detailed structure 
of V mt {7,} is unimportant W mt w 1, S mt w 0. The sys- 
tem behaves like an ideal gas and fills the volume V. At 
sufficiently low energies only the minimum of V mt {r i ] 
is explored by W mt (E — E ,N,V). The system is in a 
condensed, liquid drop (perhaps several) moving freely 
inside the empty larger volume V. 

One can guess the general form of Nln[V(E)]: Near 
the groundstate E > E it must be flat ks (N — 1) ln[Vo] + 
ln[V] because the liquid drop has some eigen-volume Vo 
in which each particle can move. With rising energy 
ln[V(i?)] rises up to the point (E trans ) where it is pos- 
sible that the drop fissions into two and Nln[V(E)} w 
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(N— 2)ln[Vo]+21n[V]. I.e. here is a first "jump" upwards 
of S mt (E) by ln{p^-}. Later further such "jumps" may 
follow. Each of these "jumps" induce a convex upwards 
bending also of the total entropy S(E) (eqEj). Each is 
connected to a bifurcation and bimodality of e s ^ E ' E ' T 
and the phenomenon of phase-separation. Of course in 
the conventional canonical picture for a large number of 
particles this is washed out and leads to the familiar 
Yang-Lee singularity of the liquid to gas phase transi- 
tion. In the microcanonical ensemble this is analogue 
to thephenomenon of multifragmentation in nuclear sys- 
tems 0,0. This, in contrast to the mathematical Y-L 
theorem, physical microscopic explanation of the liquid 
to gas phase transition shows the advantage of the mi- 
crocanonical statistics on the 3N-coordinate space as pre- 
sented here. 



B. Solid— liquid transition 

In contrast to the liquid phase, in the crystal phase a 
molecule can only move locally within its lattice cell of 
the size d 3 instead of the whole volume Vo of the conden- 
sate. I.e. in eq.lO instead we have S mt -» ln-^]^ 1 }. 

C. Summary of section Hill 

The gas-liquid transition is linked to the transition 
from uniform filling of the container volume V by the 



gas to the smaller eigen-volume of the system Vo in its 
condensed phase where the system is non-homogeneous 
(some liquid drops inside the larger empty volume V). 
First N — 1 , later at higher energies less and less degrees 
of freedom condensate into the drop. First one, then 
more and more dofs. explore the larger container vol- 
ume V leading to upwards jumps (convexities) of S mt (E) . 
The volume of the container controls how close one is 
to the critical end-point of the transition, where phase- 
separation disappears. Towards the critical end-point, 
i.e. with smaller V the jumps ln[V] — ln[Vo] become 
smaller and smaller. In the case of the solid-liquid 
transition the external volume V of the container has 
no influence on S mt {E). The entropy jumps are by 
AS""* oc lnV -lnd 3 . 

At the surface of the drops V mt > E a = min V lnt , i.e. 
the surface gives a negative contribution to S mt in eq. i|14|) 
and to S at energies E>Eq, as was similarly assumed in 
section JnJ and explicitely in eq.©. 
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